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A  HYBRID  METHOD  FOR  THE  SOLUTION  OF  SOME  MULTI -COMMODITY 


SPATIAL  EQUILIBRIUM  PROBLEMS 
by 

Jong-Shi  Pang 


ABSTRACT.  In  this  paper,  we  first  derive  a  unified  formulation  of 
the  multi-commodity  transportation  and  transshipment  spatial  equilibrium 
models  as  a  linear  complementarity  problem  with  certain  block  structure. 
We  then  show  that  a  block  successive  overrelaxation  method  is  applicable 
for  solving  the  resulting  complementarity  problem.  The  method  consists 
of  solving  a  sequence  of  subproblems  of  the  single-commodity  type.  These 
subproblems  are  to  be  solved  by  a  special-purpose  principal  pivoting 
algorithm  developed  in  an  earlier  paper.  Computational  experience  of 
solving  some  fairly  large  problems  by  the  proposed  hybrid  method  is 
presented. 


Key  Words.  Successive  overrelaxation.  Spatial  equilibrium,  Transportation, 


Transshipment,  Linear  complementarity.  Special  structure,  Computational 


experience . 


1.  INTRODUCTION 


In  tho  pas c  tow  years,  there  have  been  several  published  studies 
on  the  numerical  solution  of  various  economic  spatial  equt librium 
models  [1,  J,  5,  7,  9],  Basically,  such  models  are  concerned  with  the 
determination  of  the  regional  production,  consumption  and  inter-regional 
flows  of  several  c ommod i t i e s  so  as  to  achieve  price  oqul librium.  These 
models  have  important  practical  applications ,  see  [ -» ,  t>,  9,  II]  and 
references  therein. 

Tlte  objective  of  this  paper  is  to  propose  a  hybrid  method  tor  solving 
two  types  of  multi -conmodlty  spatial  equilibrium  models;  namely  the 
transportation  and  transshipment  models.  These  models  can  be  viewed  as 
extensions  of  their  linear  analogs  which  are  tho  well-known  linear 
transportation  and  transshipment  models.  The  method  operates  by  solving 
a  sequence  of  snbproblems  of  the  s ingle-commodl ty  type.  These  subproblems 
are  to  be  solved  by  a  spec  la l -purpose  principal  pivoting  algorithm 
described  in  an  earlier  paper  l']. 

It  Is  a  rather  well-known  fact  that  most  spatial  equilibrium  models 
can  bo  formulated  as  complementarity  problems,  see  18]  o.g.  With  the 
assumption  of  linear  supply  and  demand  functions ,  —  these  become  linear 
complementarity  problems.  Under  further  assumptions  on  the  supply  and 
demand  functions,  the  latter  complementarity  problems  can  in  principle, 
be  solved  by  a  number  of  pivoting  methods  (.most  notably,  Lemke's  almost 
complementary  pivoting  algorithm!.  However,  due  to  the  fact  that 
typically,  these  linear  complementarity  problems  are  onourmously  large, 

1/  Such  an  assumption  is  made  very  often  in  practice,  see  [0,  11]. 
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it  is  rather  impractical,  if  not  impossible  to  apply  the  pivoting  methods 
in  a  straightforward  manner. 

Among  the  recently  proposed  methods  for  solving  multi -commodity 
spatial  equilibrium  problems  the  most  promising  approach  seems  to  be 
the  Bender's  decomposition  proposed  in  [9,  10],  Nevertheless,  despite 
the  rather  encouraging  computational  results  reported  in  the  reference 
for  solving  various  practical  applications,  there  is  a  basic  drawback 
of  the  approach,  namely,  it  requires  solving  quadratic  subprograms 
(in  addition  to  solving  linear  ones)  which  could  still  be  fairly  large. 

In  essence,  the  hybrid  method  proposed  is  a  specialization  of  the 
block  successive  overrelaxation  (SOR)  Iterative  method  (described  in  [2]) 
applied  to  a  linear  complementarity  problem  with  certain  block  structure. 
As  we  shall  demonstrate,  such  a  complementarity  problem  provides  a  unified 
formulation  for  both  the  transportation  and  transshipment  spatial 


equilibrium  models. 
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rHE  rwu  spAruL  equilibrium  models 

Wo  first  dose ri bo  the  transportation  spatial  equilibrium  model.  A 

certain  number  m  of  commodities  are  to  be  shipped  from  n^  supply  locations 

to  tv,  demand  locations.  rho  problem  is  to  determine  a  set  of  flow 

variables  X  -  (X  ) .  supply  quantities  x  «  (x  demand  quantities 

«  p  i  or  1 

y  -  (y,,'.  optimal  market  supply  prices  X1  -  (\]  )  and  optimal  market 

3  1  dt  1 

2  > 

demand  prices  X'  -  (.Xg  t  3  so  that  the  equilibrium  conditions  below  are 
satisfied  for  each  a  -  1,...,n  9  -  1 . .  and  i  -  1 . m: 

(a)  Linear  regional  supply  and  demand  functions 


a  i 


m 

+  £ 
k-1 


C  V 

S0f  i  k  or  k 


■0  i 


9  i 


m 

-  £ 
k-1 


h9  l  ky9  k 
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. 


(d)  Optimum  spatial  allocation 


Ca9  1  +  Xari  '  X9i  °  ’  8  1  ■  °  ’  (C#8  i  +  X»i  '  Xpi)XaS  i"  0 


In  condition  (a)  above,  s  ,  d  ,  g  and  h  are  given  constants 

Ofi  9 1  or  i  k  9  l  K 

defining  the  supply  and  demand  functions.  In  condition  (d') ,  c  .  .  is  the 

dr  9  1 

unit  transportation  cost  from  supply  ar  to  demand  3  .  More  details  on 
the  description  of  this  model  can  be  found  in  [11]. 

1  ■> 

By  using  condition  (a)  to  eliminate  the  variables  p  and  p~,  ,  we 

or!  31 

may  formulate  the  model  as  the  linear  complementarity  problem:  find 
1  2 

vectors  x,  y,  X,  \  and  X  such  that 
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Here 


Is  the  multi-commodity  transportation  constraint  matrix,  i.e.. 
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Basically,  the  transshipment  spatial  equilibrium  model  is  defined  by 

a  set  of  equilibrium  conditions  similar  to  (a-)  -  (d)  above.  In  this  model, 

the  commodities,  Instead  of  being  shipped  directly  from  the  supply  to  the 

demand  locations,  may  go  through  some  intermediate  transshipment  points. 

If  n  is  the  total  number  of  regions  in  the  market,  then  the  problem  is  to 

find  a  set  of  flow  variables  X  -  (X  -  net  flow  quantities  y  ■  (y  . )  so 

or  0  1  cifi 

that  the  equilibrium  conditions  below  are  satisfied  for  each  a,  9  »  1 . n 

and  1  ■  1 , . . . ,m: 

(A1  Linear  regional  supply  and  demand  functions 

p«  ■  *«  -J,  l,u  V 

(B)  Conservation  of  flows 


(C)  Optimum  spatial  allocation 


ca9  i  +  Pai  '  P31 


X  .  .  >  0 

a  9  i  - 


(ca3i+Pai'p9i)Xa9i  "  ° 


Here,  we  have  assumed,  for  the  sake  of  simplicity,  that  there  is 
a  route  Joining  any  two  regions.  For  more  details  on  the  description 
of  this  model,  see  [3]. 

By  using  conditions  (A)  and  (B)  to  eliminate  the  variables  y  and  p  , 

an  ari 

we  obtain  the  linear  complementarity  problem  formulation  of  the  model: 
find  a  vector  X  such  that 
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Here  P  is  the  tnulti -commodity  node-arc  incidence  matrix^7  of  a  simple, 
complete  digraph  G  with  n  nodes  and  B  is  the  block  diagonal  matrix 


where  each  B  -  (b  )  is  m  by  m.  See  [7]  for  more  details  on  the 
ar  or  i  j 

derivation  of  the  problem  (2.3). 


2/  It  Is  obtained  by  replacing  the  ones  in  an  ordinary  node-arc  incidence 
matrix  by  identity  matrices  of  order  m.  (cf.  (2.2a)  and  (2.2b).) 
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3.  A  UNIFIED  FORMULATION 

In  the  formulation  (3.13  and  (3.33,  the  vector  X  is  arranged  according 
to  routes  between  regions;  l.e.,  conmodities  passing  through  the  same  route 
are  ordered  consecutively.  For  our  purpose,  it  would  be  more  convenient 
to  reformulate  the  problems  with  the  vector  X  partitioned  according  to 
conmodities .  For  the  transshipment  problem  (3.33,  this  can  be  easily  done. 
In  fact,  the  resulting  reformulation  has  the  form 
(3.13  Y-c  +  p:a  +  prAPX^O  ,  X^O  and  Yrx-0. 

Here  P  is  the  block  diagonal  matrix 


where  each  is  the  (single-commodity)  node-arc  incidence  matrix  of  a 
simple,  complete  digraph  with  n  nodes  and  A  is  the  block  matrix 


i 

t 


I 


A 

ran 


where  A  is  the  n  bv  n  diagonal  matrix  whose  ar-th  diagonal  entrv  is  b 

a  i  j 

The  vectors  c,  Y  and  X  are  the  corresponding  rearrangements  of  c,  Y 
and  X  respectively. 

It  is  also  easy  to  reformulate  the  problem  (3.13  in  terms  of 
commodi ties .  However,  it  is  our  contention  that  the  problem  can  be  cast 
in  the  form  (3.1)  as  well.  In  what  follows,  we  show  how  this  is  done  under 
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three  assumptions: 

(.1)  Both  matrices  S  and  D  are  nonnegative 

(2)  Each  column  of  D  has  at  most  one  nonzero  entry 

(3)  The  vector  c  is  strictly  positive. 

Observe  that  the  matrices  S  and  D  are  not  required  to  have  the 
transportation  structure  (2.2a,  b)  to  satisfy  the  first  two  conditions. 
Lemma  1 .  Suppose  that  the  assumptions  (1)  -  (3)  are  satisfied.  Then 
any  complementary  solution  to  (2.1)  must  satisfy  the  condition:  =  0. 

1  2 

Proof .  Suppose  that  (x,  y,  X,  X  ,  \~)  is  a  complementary  solution 

2 

to  (2.1)  wi th  >  0  for  some  component  i.  Then  by  complementarity,  it 

2 

follows  that  V  =  0.  Let 

J  -  Jj:  DtJ  t  0)  . 


Then  for  each  j  € J,  we  have 


Yj  '  Cj  +  (sTx1)j  ‘  DijXi  =  Cj  +  >  0 


which  implies  that  Xj  =  0.  Hence, 


which  is  a  contradiction. 


Q.E.D. 


Referring  to  the  transportation  spatial  model.  Lemma  1  says  that  if 
the  transportation  costs  are  strictly  positive,  then  no  excess  demand  is 
possible  at  equilibrium. 


10 


Proposition  2.  Under  the  same  assumptions,  problem  (2.1)  is  equivalent 
to 

(3.3)  u=s  +Gx  -X>0  x  >  0 

=  =5 

Y  =  (c  +  DTd)  +  DT  H  D  X  +  STX  >  0  X>0 


x  -  SX 


>0  X  >  o 


Tit 

u  x  =  Y  X  =  m.  X  =  0  . 


in  the  sense  that  a  complementary  solution  to  one  problem  will  always 
give  rise  to  one  such  solution  of  the  other. 

1  2 

Proof .  Let  (x,  y,  X,  X  ,  X  )  be  a  complementary  solution  to  (2.1). 
Then  by  Lemma  1 ,  we  have  y  =  DX.  Hence 


Y  =  c  +  sV  -  DT(v  -  d  -  Hy) 


or  equivalently, 

(3.4)  Y  +  DTv  =  (c  +  DTd)  +  DTHDX  +  STx’  . 

We  claim  that  (x,  X,  xS  is  a  complementary  solution  to  (3.3).  It  suffices 
to  show 


XT(c  +  DXd  +  DTHDX  +  sV)  =  0. 


According  to  (3.4),  the  left-hand  term  is  equal  to 


XT(Y  +  DXv)  -  (DX)Xv  -  yXv  -  0 


Conversely,  suppose  that  (x,  X,  X)  i9  a  complementary  solution 


to  (3.3) .  Define 


n 


v  -  \  .  v  -  n\  an, I  \'  -  v>'  •  MV' 


»  t»  I  ha  llaflat  tVa  part  of  tha  vaotot  ?.  1,a  I'lolw  that  ^\.  V.  \, 

\  ,  \-'  t«  .1  .  .Mliplamallt  .11  Y  Bolllt  toil  tO  V*  ■*'  r,<  «llOV  (hat  f  t  1  a 

faaathla,  wa  hava 


'  , 

*t  ♦  Uy  ♦  \ '  -  \,i  *  uv'  o 


vtlaia  ?  t«  tha  \V»UU-(1  \'AVt  of  tha  Va.tOt  9  ,  MoVanVai  .  It  t  «  ol'VtOI," 
that  tha  follovtus  Itol. In 


Y  >  v  ' 

.■  ♦  sV  It  \' 


-  v>-  •  p\n  >  h'ii nx  •  s ’ \ '  i'\.i  •  Mvt 


If  V<-  *  S  \  I'  \  '  ^  •  i'  Cot  »oma  .owi'onant  |,  (lieu  .olimtn  I  of  I'  wuat 
.ontatw  *«i'(  Iy  .'tin  poatttva  ant  iv  ,  aav  p  -  i' .  Cm  <  havmot a  ,  va  mnat 
have  \  .1  1  IIV  '  ^  -  0  I «tp  I V  t  \Y)i  \j  -  |'  Hall,' a 

"  •  sV  "v\  -  *  sV\  >',,'1 

vlil.li  t  «  tmppaa  tl'lts,  I'lusi  afot  t»  t  |i«s  vactot  i\.  v.  \.  v'  1®  (rntalMa 

to  v I'.  IV  v.'tiH'Uta  tha  pivof,  It  rama tn«  to  aliov  (list 


r  j  !•  v  1  r  ' 

Y  v»l  *  'tv'  -  \  VO  *  S  v  |i  V  '  -  0  . 


Snppoae  that  v«l  «  Ifcr'  '  t'.  It  than  follovn  that  Vj  -  0  .  ltv  iialiig  an 
ai)iiniiatit  •  initial  to  tha  proof  .'f  I. amnia  1.  w  wav  raa.ltlv  .ta.lo.a  that 
V(  -  0.  t'onaa.jnant  Iv  .  It  hoi, la  that  v\.l  *  Mv '  *  -  t' .  Filial  Vv,  It  follows 
t  tv  mu  v'  "  that 

v 1  v «  *  ;•  v '  o'\''  -  \'u«'  •  h'.i'  *  i'rui'\  *  s’\'|  fi'\'\,i  *  i^'* 
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Remark.  TT\e  second  part  of  the  proof  remains  valid  If  c  is  merely 
nomxegative. 

Ihe  last  proposition  establishes  a  first-stage  transformation  of 
the  problem  (2.1).  To  complete  the  process,  we  prove 

Lemma  3 .  If  the  matrix  S  is  nonnegative,  then  any  complementary  solution 

to  (3.3)  must  satisfy:  uT(SX)  -  -  0. 


Proof.  In  fact,  if  u^  >  0,  then  x.  »  0  by  complementarity.  Hence 

0  i  -  -  (SX)t  v  0 

from  which  the  desired  equalities  follow  readily.  Q.E.D. 


Proposition  4.  Suppose  that  assumptions  (1)  -  (31  are  satisfied.  Then 
problem  (3.31  (and  thus  (2.111  is  equivalent  to 

(3.61  y  -  (c  +  DTd  +  STsl  +  (DTHD  +  STGS1X  +  STG^  >0  ,  X  >  0 

X-  s  +GSX  +  G(i  >0  ,  u  >  0 

T  T 

n>\n*o 

in  the  sense  specified  in  Proposition  2.  The  latter  linear  complementarity 
problem  can  be  written  in  the  form 

T  T  T 

w  -  r  +  P  a  +  P  BPz  £  0  ,  z  >  0  wr.  -  0 
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where 


Prooj;.  Suppose  chat  (x,  X,  \)  is  a  complementary  solution  to  (3.3). 

We  claim  that  (X,  p)  la  a  complementary  solution  to  (3.t>).  in  fact, 
we  have 

a  +  CSX  +  C  n  *,s+Gx«u  +  \'>.o 

and 

(c  +  DTvl  +  $rs)  +  (DTHD  +  STGS)X  +  STo'u 
-  (c  +  DTd)  +  DTHa\  +  ST\  +  sru  >  0 

m 

establishing  the  feasibility.  The  desired  complementarity  relationships 
follow  from  Lemma  J. 

Conversely,  let  (X,  pi  be  a  complementary  solution  to  (3.61 .  Define 
x  ■  u  +  SX  and  \  ■  s  +  Gx  . 

Obviously  (x,  X,  \)  solves  (3.3)  .  Q.E.D. 

The  last  proposition  shows  that  if  the  transportation  costs  are 
all  positive,  then  the  transportation  spatial  equilibrium  problem  (2.1) 
can  be  cast  as  a  linear  complementarity  problem  of  the  same  type  as  the 


’’1 
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one  obtained  from  the  transshipment  model.  In  particular,  by  partitioning 
the  flow  vector  (XI  and  the  excess  supply  vector  in')  according  to 
commodities,  we  may  rewrite  the  linear  complementarity  problem  (3.01 


in  the  form  (3.1)  and  (3.2)  with  each  and  A  given  by 


(3.7)  Pt  - 


Di  ° 


Sl  1 


I  AU  ° 


o  a* 


i  ■> 

where  A  and  A  are,  respectively,  n,,  by  n,  and  n  by  n  diagonal  matrices  with 

lj  LJ  i  i  \  \ 


(AIj)or«  "  ha  ij  and  ^>9  3  "  *B  U  ;  and  where  (. 


is  the 


s ingle -commod i ty ,  transportation  constraint  matrix,  i.e.,  (2.2a)  and 
(2.2b)  with  m  ■  1 . 

It  is  interesting  to  note  that  in  both  cases,  the  matrices  A 
are  diagonal  and  the  P^  are  identical.  In  fact,  each  P  is  the  node-arc 
incidence  matrix  of  the  graph  (or  network)  underlying  the  model.  This 
is  clear  in  the  transshipment  case.  In  the  other  case,  the  P^  (cf. 

(3.7))  can  be  considered  as  the  node-arc  incidence  of  a  transportation 
bipartite  graph  with  a  self-loop  connecting  each  supply  node.  Such  a 
loop  represents  the  excess  supply  at  that  node. 

To  summarise,  we  conclude  that  both  the  multi -conmodlty 
transportation  and  transshipment  spatial  models  can  be  formulated  as 
a  linear  complementarity  problem  of  the  form  (note  the  change  of  notation) 


(3.8)  w  -  r  +  PTs  ♦  PTAPz  >  0 


z  >  0 


wTz  »  0 


where  the  matrix  P  is  block  diagonal  with  the  number  of  blocks  equal  to 
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the  number  of  comnodltles .  We  point  out  that  the  matrix  pTAP  has 
the  block  structure 


(3.91 


/X  Au p. 


T 

Pm  A  P, 
m  ml  1 


-  -?]  A  P  \ 

!  lm  m  i 


-  -P  A  P 
m  mm  m 
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4.  THE  HYBRID  METHOD 

There  are  cvo  basic  reasons  motivating  the  proposal  of  a  hybrid 
method  for  solving  the  two  spatial  models  described  earlier.  The  first 
is  the  large  size  of  the  resulting  linear  complementarity  problem  which 
is  typical  in  practice.  For  instance,  a  transshipment  problem  with 
1 0-commod l t ies  and  50-regions  would  give  rise  to  a  24,500  by  24,500  matrix 
in  the  complementarity  problem.  The  second  reason  is  due  to  the  recent 
development  of  an  efficient  special-purpose  algorithm  for  solving  the 
single-commodity  models-  (see  [7]>.  It  is  hoped  that  this  algorithm 
would  provide  an  effective  subroutine  for  solving  the  multi-commodity 
problems . 

To  state  the  hybrid  method,  we  refer  to  the  unified  formulation 
of  the  two  models,  i.e.,  the  linear  complementarity  problem  (3.8")  and 
(3.9).  Recall  also  chat  the  matrix  P  is  block  diagonal  with  diagonal 
blocks  . 

Algorithm 

Ste£__0.  Let  z  *  (*^,...,2  )  be  an  arbitrary  nonnegative  vector  and 
let  I  be  a  scalar  with  0  <  i  ^  2.  Put  k  ■  0  . 

•  For  i  ■  1,...,m,  solve  the  linear  complementarity  problem 


3/  Although  the  reference  treats  only  the  single-commodity  transshipment 
model,  due  to  the  fact  that  the  transportation  model  gives  rise  to 
exactly  the  same  type  of  linear  complementarity  problem  it  is  not 
difficult  to  extend  the  algorithm  described  in  the  reference  to  handle 
the  transportation  model.  See  section  5  for  computational  results  on 
this  extension. 


u 
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(4.1)  m  rt  +p[s*  +  p[A.iiPizi  >0  ,  2t>0  and  w^  =  0 

4/  lC  lr4-1  I#  I.  i  1 

with-  s  ■  s  +  E  A.  ,P,Z,  +  ^  A.  ,P.Z,  »  *et  z.  *  be  a  solution 

J  <1  J  J  j  >i  J  J  1 


and  set 


k+1  k  k+1 .  k+g  k  . 
*1  "  Zi  +  “1  (21  ‘  *1  > 


where 


k+1 


t  -  “  k+1  „  < 

max  (iu  :tu  <  u>  ,  z  >  0  J 


Step  2.  Let 


C(i.J):  <zi+1> 


>0  ,  or  (w^+1)j  -  O^  +  pJ^ 


m 

•  I 
•L-1 


Anpt*r»j  <o> 


If  max  I(w^+')«  I  S  •  where  «  is  a  given  positive  tolerance,  stop. 

(1,J)6J  1  J 

An  approximate  solution  is  at  hand.  Otherwise  return  to  Step  1  with  k 
replaced  by  k+1 . 

In  essence,  the  algorithm  is  a  specialization  of  the  modified  block 
SOR  method  for  solving  a  general  linear  complementarity  problem  with 
block  structure  (see  [2]).  Concerning  its  convergence,  we  state  the 
following  principal  result. 


Theorem  5.  Suppose  that  the  matrix  A  is  symmetric  positive  definite. 
Suppose  also  that  for  i  -  1 . .  the  system 


(4.2)  0  <  yt  +  0  ,  ?iyi  -  0 


T 

Vi 


<  0 


4/ 


The  vectors  r  and  s  are  partitioned  as  r  -(r*, 
in  accordance  with  the  matrix  P. 


,r£)T  and  s  »  (s*. 


.v 

m 
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is  inconsistent.  Then  the  sequence  £zk  -  (z^ . z^)}  of  vectors 

generated  by  the  algorithm  has  a  limit  point.  Moreover,  any  such  limit 
point  solves  the  linear  complementarity  problem  (3.8). 

It  should  be  pointed  out  that  even  though  the  matrix  A  is  positive 
definite,  the  matrix  PTAP  (and  each  of  its  diagonal  submatrices 
p^A^P^)  may  not  be  so.  Hence,  the  theorem  is  not  a  direct  consequence 
of  the  convergence  result  established  in  the  reference.  Nevertheless, 
it  can  be  proved  in  very  much  the  same  way.  In  what  follows,  we  sketch 
the  proof. 

First  of  all,  the  inconsistency  of  the  systems  (4,2)  and  the  positive 

definiteness  of  the  matrix  A  imply  the  feasibility  and  hence  the  solvability 

k+i 

of  each  of  Che  subproblems  (4.1).  Thus  each  vector  zt  2  is  well-defined. 

Next,  by  using  the  function  f(z)  ■  (r  +  P^s)^z  +  $(Pz)^APz  to  monitor  the 

progress  of  the  algorithm,  it  is  not  difficult  to  show  that  the  sequence 

{f(z  )}  is  monoConically  decreasing.  From  this  it  follows  chat  the 

sequence  {z^}  must  be  bounded  and  therefore  has  a  limit  point.  Finally, 

that  any  such  limit  point  solves  the  linear  complementarity  problem  (3.3) 

can  be  proved  by  using  the  positive  definiteness  of  the  matrix  A. 

T 

Without  the  positive  definiteness  of  the  matrix  P  AP ,  there  is  in 

k 

general  no  guarantee  that  the  entire  sequence  {z  }  will  in  fact  converge. 
However,  the  next  result  shows  that  the  sequence  (Pz  ]  will. 

Proposition  6.  Under  the  assumptions  of  Theorem  5,  the  sequence  {Pz  } 
converges  to  Pz*  where  z *  is  a  limit  point  of  {z  }  . 


proof.  Observe  that  any  limit  point  of  the  sequence  {Pz  }  must  be  of 
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the  fora  Pz*  where  z*  is  a  limit  point  of  the  sequence  fzk]  .  The 

fact  that  {Pz  }  has  at  least  one  limit  point  follows  from  its 

k 

boundedness  (which  is  implied  by  the  boundedness  of  {z  }).  Suppose  now 
that  {Pz^}  has  another  limit  point  given  by  Py*  where  y*  is  some  limit 
point  of  {z  }  .  By  Theorem  5,  both  y*  and  z*  solve  the  linear 
complementarity  problem  (3.8).  By  the  positive  definiteness  of  A,  it 
can  be  proved  easily  that  Py*  -  Pz*.  Thus,  the  sequence  {pzk}  has 
a  unique  limit  point  and  therefore  converges.  Q.E.D. 


Specializing  Theorem  5  to  the  two  spatial  models,  we  obtain 

Theorem  7.  Suppose  that  each  of  the  matrices  G  .  H.  and  B  (cf.  (2.2c) 

at  s  at 

and  (2.4))  is  symmetric  positive  definite.  Suppose  also  that  the  vector  c  of 
transportation  costs  is  positive  (in  both  models).  Then  the  assumptions 
of  Theorem  5  are  satisfied.  In  particular,  the  sequence  of  vectors 
generated  by  the  algorithm  will  converge,  in  the  sense  specified  in  Theorem  5 
and  Proposition  6,  to  a  solution  of  the  two  models. 

The  proof  of  Theorem  7  is  easy  and  thus  omitted.  Observe  that  in 
the  specialization,  each  subproblem  (4.1)  is  precisely  one  of  the 
single-commodity  type  to  which  the  algorithm  described  in  [7]  can  be 
applied. 
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5.  COMPUTATIONAL  EXPERIENCE 

In  this  section,  we  report  our  computational  experience— ^  with 
the  hybrid  method  for  solving  some  randomly  generated  spatial 
transportation  and  transshipment  models.  For  completeness,  we  first 
present  a  brief  summary  of  the  computational  results  (see  Tables  1 
and  2)  with  the  special-purpose  principal  pivoting  algorithm  (described 
in  17])  for  solving  the  s ingle -comnodity  problems.  Me  explain  how  the 
data  were  generated.  For  the  transportation  model  (cf.  the  formulation 
(3.6)),  the  transportation  costs  were  randomly  generated  in  the  interval 
(0,10),  the  components  of  the  vectors  s  and  d  in  (-50,0)  and  (-  10,0) 
respectively  and  the  diagonal  entries  of  the  matrices  G  and  H  in  (0,5). 

The  data  for  the  transshipment  model  were  generated  in  the  same  way  as 
in  [7], 

The  third  column  (basic  arcs)  in  Tables  1  and  2  refers  to  the 
number  of  flow  variables  that  are  positive  in  Che  final  solution. 

Notice  that  this  number  is  bounded  above  by  the  total  number  of  regions  less 
one.  The  fourth  column  (excess  supplies)  in  Table  1  refers  to  the 
number  of  supply  regions  where  the  actual  supply  quantity  exceeds  the 
effective  supply  shipped. 

As  with  all  relaxation  methods,  the  efficiency  of  the  hybrid 
method  proposed  is  critically  dependent  on  the  choice  of  the  parameter  u>  . 

In  the  experimentation,  we  have  tried  several  values.  In  Tables  3  and  4 
below,  we  report  a  summary  of  the  results  with  the  use  of  the  particular 

5/  All  the  computations  were  performed  on  a  DEC-20  computer  and  in  double 
precision.  The  computer  codes  were  written  in  FORTRAN.  All  the 
timings  reported  are  exclusive  of  input  and  output. 
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(#suppl les , 
^demands) 

(Fares 

(Fbasic  arcs 

(Fexcess 

supplies 

(Fpivots 

CPU  time 
Total 

(in  sec.) 
Per  pivot 

(80,801 

6,400 

100 

57 

794 

298.983 

0.389 

(80,40) 

3,200 

55 

65 

399 

77.307 

0.194 

(60,120) 

7,200 

138 

38 

933 

406. 102 

0.435 

(60,100) 

6,000 

129 

27 

707 

250.170 

0.354 

(40,120) 

4,800 

148 

4 

603 

180.451 

0.299 

(20,100) 

2,000 

96 

1 

192 

26.388 

0.  137 

Table  1  . 

Single-conmodi ty 

transportation  spatial 

.  model 

(Fregions 

(Fares 

(Fbasic  arcs 

(Fpivots 

CPU  time 
Total 

(in  sec.) 
per  pivot 

120 

14,280 

71 

377 

158.397 

0.420 

105 

10,920 

62 

142 

45.527 

0.321 

90 

8,010 

59 

31 1 

76.272 

0.249 

75 

4,150 

43 

175 

29.232 

0.167 

30 

870 

17 

35 

.945 

0.027 

Table  2.  Single -commodity  transshipment  spatial  model 


x  jiving  the  fastest  convergence.  The  data  for  the  problems  solved 
were  generated  in  the  same  w*v  as  in  the  single-commodity  case.  The 
tolerance  c  for  the  termination  criteria  (cf.  Step  2  in  the  algorithm) 
was  chosen  to  be  10  '  .  As  the  task  of  testing  for  termination  in 
each  iteration  is  rather  time-consuming  (because  of  the  large  number 
of  variables  in  the  problem) ,  we  chose  not  to  perform  such  test  at  every 
iteration,  but  rather  at  every  5  iterations.  The  entries  in  the  two 
tables  below  should  be  se If -explanatory . 


’.‘supplies 

‘.‘demands  **c 

ommodi ties 

20 

20 

3 

15 

15 

9 

15 

15 

6 

10 

10 

9 

15 

15 

3 

10 

10  6 

Table  3.  Multi -Cotnmodi  ty 

‘‘regions 

•^commodities 

Ol 

40 

<> 

1.7 

20 

10 

1.7 

20 

6 

1.7 

15 

10 

1.7 

15 

4 

1.7 

10 

10 

1.75 

10 

s 

1.75 

total 


01 

‘^iterations 

CPU  time 

1 .4 

30 

264. S71 

1  .4 

35 

309.769 

1  .4 

55 

371 .865 

1.4 

35 

98. 998 

1 .2 

15 

56 .4S7 

1  .2 

30 

02.253 

transportation  spatial  model 
^iterations  Total  CPU  time  (sec.) 


140 

592.419 

125 

326.102 

120 

150.913 

60 

56.449 

45 

18.355 

85 

29.799 

35 

3.886 

Table  4. 


Multi-commodity  transshipment  spatial  model 


<3.  CONCLUSION 


In  this  paper,  we  have  shown  how  a  hybrid  iterative  method  can  be 
used  to  solve  the  multi-consnodity  transportation  and  transshipment 
spatial  equilibrium  models.  At  this  stage,  we  are  unable  to  compare  the 
numerical  performance  of  the  proposed  method  with  other  methods  that  are 
applicable.  To  the  best  of  our  knowledge,  there  is  no  published  numerical 
results  for  solving  problems  as  large  as  the  ones  solved  in  this  paper. 
(The  largest  problem  solved  by  the  Bender's  decomposition  approach  has 
only  eight  supply  and  demand  regions  and  four  commodities .  See  [9,10].) 

It  would  be  interesting  to  see  how  these  problems  can  actually  be  solved 
by  other  methods,  especially  the  Bender's  approach. 
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